Modulation theory of quantum tunneling into a Calogero-Sutherland fluid. 
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Quantum hydrodynamics of interacting electrons with a parabolic single particle spectrum is stud- 
ied using the Calogero-Sutherland model. The effective action and modulation equations, describing 
evolution of periodic excitations in the fluid, are derived. Applications to the problem of a single 
electron tunneling into the FQHE edge state are discussed. 



It is known since 60's[l[ that Landau Fermi liquid the- 
ory is not applicable to one dimensional metals. The 
thermodynamic of one dimensional electronic system is 
explained by Tomonaga-Luttinger (TL) modeldH. This 
model is a simplified version of a quantum hydrodynam- 
ics that takes into account only a single mode of sound 
excitation (phonons) of an electronic fluid. This ap- 
proximation is equivalent to replacing a true electronic 
single-particle spectrum by the linear spectrum of Dirac's 
fermions. Though the TL model accounts for thermody- 
namics it fails for effects that involve breaking of the 
particle-hole symmetry. In other words only the diago- 
nal part of kinetic coefficient matrixQ can be found using 
this approach, while all non-diagonal terms vanish in this 
approximation. 

Among phenomena missed by the TL approximation 
are thermopower, photovoltaic effect and Coulomb drag 
with a small momentum transfer. To solve any of this 
problems one needs to use quantum hydrodynamic with- 
out making one mode approximation. 

The problem of single particle spectrum curvature and 
electron-electron interaction was recently addressed in 
Ref.®, who proposed to mimics a nonlinear spectrum 
by two types of electrons with different Dirac's spectra. 
This approach agrees with perturbative expansion of ex- 
act results for Calogero- Sutherland (CS) model 0. 

Here, we would follow a different route of "exact" 
bosonization, originally developed in Ref.® for the 
sake of matrix models. As expected®, a curvature of 
fcrmionic spectrum leads to the cubic in density terms in 
a bosonic Hamiltonian. Though formally a bosonic de- 
scription is achieved, the resulting theory is non linear. 
This renders this method extremely difficult for appli- 
cations and no new results based on it are reported, as 
yet. 

To gain an understanding of one dimensional hydro- 
dynamics we study the special case of electrons with sin- 
gular 1/r 2 interaction. This problem is known as the 
Calogero- Sutherland model (CS)® and has numerous 
applications in condensed matter and nuclear physics (lfjj|. 
Being exactly solvable 11 , CS model posses an infinite 



number of integrals of motion. That render its hydrody- 
namics special and easier to analyze. 

The connection between the hydrodynamic excitations 
of the CS fluid and correlation functions of various opera- 



tors is not straightforward. The simplest known example 
is the soliton, a solution of hydrodynamical equation of 
motion propagating without changing its shape. In orig- 
inal description it corresponds to the many body exci- 
tation, called anyon. Anyons obey a special "fractional" 
statistics, which is neither Fermi nor Bosc like. The cor- 
relation function of anyons is similar to the one of free 
fermions, with the values of Fermi velocity and momen- 
tum renormalized by interaction 13 1. 

In order to find other correlation function (such as elec- 
trons' Green function) one needs to solve the system of 
saddle point equations, that are the continuity and the 
Euler equation for CS fluid. This is a very difficult prob- 
lem. It is considerably simplified if the solution is a peri- 
odic wave with slowly changing wave parameters. In this 
case, the original hydrodynamic equations may be re- 
placed by Whitham's modulation equations 14[. To put 
it simple, the TL model describes particle-hole excita- 
tions of fermionic problem as sound modes (vibrations) 
of one dimensional harmonic string. The modulation the- 
ory describes particle-hole excitations as vibration of an 
anharmonic string. 

The modulation technique is well developed for various 
models of dissipationless hydrodynamics. In particular 
for the Benjamin-Ono equation ]WL flt3 |. that describes 
dynamics of internal waves in stratified fluids of great 
depth. As it was stressed recently [17}, this two problems 
are mathematically related. Remarkably, the classical 
BO equation appears in FQHE[18|, governing the evo- 
lution of a semiclassical wave packets containing a large 
number of fermions. It is therefore not surprising that 
some formal aspects of the derivations of modulation the- 
ory for the BO and the CS models are similar. The final 
results of these two models are nevertheless different. 

The structure of this work is a following: starting with 
the microscopic description of the CS model we pass into 
its hydrodynamics (the technical details are given in the 
Appendix). Next, we consider periodic density excitation 
of the CS model and derive its evolution equations. Fi- 
nally, we apply this machinery to the problem of a single 
electron tunneling into the CS liquid. 

The microscopic Hamiltonian of CS model is given by 
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Here A is a strength of particle interaction, m is an elec- 
tron mass, Xi is a coordinate of the particle on a circle 
with the perimeter L. Now on we use the convention 
h = 1, m = 1. 

On large scales and long time CS model is described 
by hydrodynamical Hamiltonian 
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(the details of derivation that follows Ref. 0, are pre- 
sented in Appendix). The whole specific of CS model is 
incorporated in potential energy term 
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where the Hilbert transform is defined as 
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The conventional TL model can be obtained by expan- 
sion of the Hamiltonian © up to second order in fluctua- 
tions (due to conservation laws the velocity v and density 
operators bp are of the same order) and neglecting the 
high gradient terms. 

It is convenient to reformulate the problem in La- 
grangian description, defining the action 
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This action reaches its minimum provided saddle point 
equations are satisfied: the continuity 
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and Eulcr equation 
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Here the enthalpy 
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The dependence on the interaction constant in the limit 
of strong interaction (A 3> 1) had been eliminated by the 
rescaling (x = x/X,p = p\,i , = t/X). Now on we use the 
rescaled coordinates, omitting the tilde. 

To derive modulation equation it is convenient to in- 
tegrate the velocity field out 



L 



dx 



2p 



H , 



6 



P + „FA 



8 p 



(9) 



This Lagrangian has one-periodic density wave solutions 
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where 6 — kx — Lot. The dispersion relation between 
amplitude of the wave 
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To proceed further we replace the Polychronakos solution 
with a modulated one 



9{x, t) = k{x, t)x — lu(x, t)t 



(13) 



and allow po to depend on coordinate and time as well. 
Strictly speaking this is no longer the minimum of the 
action ([5]) . However under the condition the modulation 
technique works, this solution minimizes the action "in 
average" . We define Lagrangian averaged over one period 
of oscillations 
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After some straightforward, though lengthy calculation, 
one finds a Lagrangian of anharmonic string 
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Here 7 = — d x 1 dtPo- Applying the least action condition 
to the averaged Lagrangian 
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one obtains modulation eqs. for CS model 
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These eqs. govern the dynamics of the wave parameters 
in the density excitations propagating through CS fluid. 

Next, we apply this theory to study an evolution of 
a distortion caused by adding one electron to the fluid. 
This problem is a prototype for a quantum tunneling from 
the normal metal into the edge state of FQHE. 

An added electron causes a density fluctuation that 
splits into two chiral parts, moving in the opposite direc- 
tions. Using Riemann invariants [2 1| 
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one approximates egs. (|6I8|) by 



vector given by 



ut + uu x = 
u t + uu x = 0. 
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The tunneling at point x = corresponds to the initial 
conditions Uo(£) = > ( w here e has a scale of an elec- 
tron wave length). To study right chiral sector we pass 
into a reference frame moving with a sound velocity to 
the right (£ = x — npt). Solving Hopf cq. (f2"0"]) by Godo- 
graph method we find an implicit solution 



u = wq(£ — ut) 
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Solving cubic equation ([2"Tj) , we find an explicit solution 
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where £_ = 3(et/4) 1 / 3 , ^ + = t/e are tailing and leading 
edge coordinates. Inside the interval (£_ < £ < £+) so- 
lution of eq. pT) is multivalued and we denote three dif- 
ferent branches by f± > fi > fa. The multi- validity of 
solution reflects a wave breaking phenomena. The sin- 
gle value solution inside the interval (£_ < £ < £+) is 
restored by keeping the second-order spatial derivatives 
in eq. ijHJ . An elegant way of dealing with second-order 
gradient terms is to use the modulation technique de- 
veloped above. Let us assume that for (£_ < £ < £ + ) 
density is an oscillating function satisfying modulation 
cqs. (|17[) with a proper boundary condition 22 1. 

At the boundaries, the particle density found from 
Hopf equation, should match the density, averaged over 
the period of fluctuations, inside the oscillating interval 
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In addition, the amplitude of the oscillation vanished 
at the trailing edge 

A = 0, C = £- ■ (23) 

and the wave vector vanishes at the leading edge 

k = 0, £ = £+ . (24) 

It is convenient to define phase velocity c = uj/k and 
hydrodynamic velocity V = j/po- In a new variables 
eq.(|23[) can be rewritten as 
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Solving eas.([T7)). we find that electron tunneling into CS 
liquid excites density wave with phase velocity and wave 
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As we see, the excitation that follow electron tunneling 
are quite different from the anyon tunneling [13]. The 
later resulted in the creation of a single soliton that was 
propagating with a constant velocity through the fluid. 
Electron tunneling causes a spreading density evolution 
that consist many picks and linearly increasing phase ve- 
locity. 

The large number of oscillations experienced by the 
fluid density and the smooth dependence of wave pa- 
rameters a posteriori justifies the validity of modulation 
technique. One can check, using eq. (f2"3"]) . that solution 
eq. (f26]) respects the particle conservation and accounts 
for the half electrons moving to the right (solution for 
another half of electron moving to the left differs by sign 
change of a sound velocity). 

The periodic density wave that develops after the tun- 
neling can be viewed as a superposition of anyons, with 
individual anyone corresponding to the picks of the den- 
sity oscillations. Therefore eq. pB"]) describe a decay of 
an electron into a large number of quasiparticles of CS 
model. This process may be detected by measurement 
the time dependence of a current that follows the tun- 
neling event. Unlike the fractional charge measurements 
in shot noise experiments 23|, 24 1 , the tunneling discussed 
above is from normal metal to the FQHE state. There- 
fore the fractionalization of the elementary charge does 
not show up in a low frequency shot noise, but in the 
finite frequency noise of a current pulse. 

In this work we developed modulation theory for the 
CS model. We applied this theory to the problem of a 
single-electron tunneling. We found an evolution of a 
current pulse that followed the tunneling event. 

This work was motivated by discussion with I. 
Gruzberg and P.B. Wiegmann, from whom I learned var- 
ious mathematical ideas used in this paper. I benefited 
from discussion with A. Kamenev, D. Maslov, A. Mir- 
lin, M.Stepanov and M.Stone. I acknowledge the Memo- 
rial University of Newfoundland, where large part of this 
work was done, for the hospitality. My research was sup- 
ported by NSF-DMR-0308377. 



APPENDIX:DERIVATION OF 
HYDRO DYNAMICAL THEORY 

Consider a ring geometry. Electron's coordinate is rep- 
resented by a complex variable z n = Le % ® n where 9 is an 
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angle along the circle of radius L. In this variables the 
Hamiltonian eq.JT]) is given by 
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The ground state of the Hamiltonian wave function of 
(|2T|) can be found exactly 
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is Vandermonde determinant. Excited states are given 
by 



*k = *n J K . 
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where Jack polynomials J K are parameterized by parti- 
tion k. The problem had been reduced to the properties 
of the new bosonic Hamiltonian 
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that acts in the Hilbert space of symmetric wave func- 
tions. It is given by 
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where Di — Zidi. Aiming for a second quantization one 
defines so called collective variables 
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In terms of collective variables the bosonic Hamiltonian 
can be rewritten asEl 
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So far the transformation have been exact. Passing to 
the hydrodynamic limit (N — > oo, oo 2N/L — * p 
one arrives to eq.@; a; is a coordinate along the circle 



(x = the linear density p(x) = ^p(0). The modes 

of velocity operator are defined as 
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It is easy to see that these definitions are consistent with 
a standard commutation [25| 

[v(x), p(y)] = -iS'(x - y) . (36) 
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